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Fig. 25. Plasmodium {Cancer JPagwrus). 

Fig. 26. Corpuscles of Aster acanthion vulgare, freshly drawn. 

Figs. 27 — 28. Union of a group of corpuscles of AsteracantMon vulgare. 

Fig. 29. Portion of a plasmodium produced by tKe union of tlie finely granular 
corpuscles of JEchinus sphmra, showing distinct endoplasm containing 
the coarsely granular and the coloured corpuscles, and ectoplasm sending 
out filamentous pseudopodia, which unite with those of free corpuscles. 

Figs. 30 — 33, Fhoner gates vorax, from " Zeitsch. f. Wiss. Zool." Bd. XXX. 1878. 
Taf. 11, figs. 54-57. 
All the figures drawn with Yerick, Oc, 2, Obj. 7. 



I V, *' On tlie Analytical Expressions which give the History of 
a Fluid Planet of Small Viscosity, attended by a Single 
Satellite/' By G. H. Darwin, F.R.S. Received March 6, 
1880. 

In a series of papers read from time to time during the past two 
years before the R-oyal Society, I have investigated the theory of the 
tides raised in a rotating viscous spheroid, or planet, by an attendant 
satellite, and have also considered the secular changes in the rotation 
of the planet, and in the revolution of the satellite. Those investi- 
gations were intended to be especially applicable to the case of the 
earth and moon, but the friction of the solar tides was found to be a 
factor of importance, so that in a large part of those papers it became 
necessary to conceive the planet as attended by two satellites. 

The differential equations which gave the secular changes in the 
system were rendered very complex by the introduction of solar 
disturbance, and I was unable to integrate them analytically ; the 
equations were accordingly treated by a method of numerical quadra- 
tures, in which all the data were taken from the earth, moon, and 
sun. This numerical treatment did not permit an insight into all the 
various effects which might result from frictional tides, and an analy- 
tical solution, applicable to any planet and satellite, is desirable. 

In the present paper such an analytical solution is found, and is 
interpreted graphically. But the problem is considered from, a point 
of view which is at once more special and more general than that of 
the previous papers. 

The point of view is more general in that the planet may here be 
conceived to have any density and mass whatever, and to be rotating 
with any angular velocity, provided that the ellipticity of figure is not 
large, and that the satellite may have any mass, and may be revolving 
about its planet, either consentaneously with or adversely to the plane- 
tary rotation. On the other hand, the problem here considered is 
more special in that the planet is supposed to be a spheroid of fluid of 
small viscosity ; that the obliquity of the planet's equator, the inclina- 
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tion and tlie eccentricity of satellite's orbit to the plane of reference are 
treated as being small, and, lastly, it is supposed that the planet is 
only attended by a single satelHte. 

The satellite itself is treated as an attractive particle, and the planet 
is supposed to be homogeneous. 

The notation adopted is made to agree as far as possible with that of 
a previous paper, in which the subject was treated from a similarly 
general point of view, but where it was supposed that the equator and 
orbit were co-planar, and the orbit necessarily circular.* 

The motion of the system is referred to the invariable plane, that is, 
to the plane of maximum moment of momentum. 

The following is the notation adopted : — 

For the planet : — ■ 

If = mass ; a = mean radius ; g — mean pure gravity ; C = moment 
of inertia (neglecting ellipticity of figure) ; n = angular velocity of 
rotation ; i = obliquity of equator to invariable plane, considered as 
small ; g = f ^/c&. 

For the satellite : — 

m = mass ; c = mean distance ; O — mean motion ; e = eccentricity 
of orbit, considered as small ; j •= inclination of orbit, considered as 
small ; t =£ ^m/o^, where m is measured in the astronomical unit. 

For both together: — ■ 

1/ = M/m, the ratio of the masses; s = f[('^^'/^)^'(l + ^)]^; ^ ^== "tlie 
resultant moment of momentum of the whole system ; E = the whole 
energy, both kinetic and potential, of the system. 

By a proper choice of the units of length, mass, and time, the nota- 
tion may be considerably simplified. 

Let the unit of length be such that M~\- m, when measured in the 
astronomical unit, may be equal to unity. 

Let the unit of time be such that s or f[<^'^/^)'^(l-i-^)]^ may be 
unity. 

Let the unit of mass be such that C, the planet's moment of 
inertia, may be unity. 

Then we have 

Q3c3 = if-|-m=:l ....... (1). 

ISTow, if we put for g its value M/a^, and for v its value if/m, we 
have 



2 



5 



9- 



'aM w 



m M 



_______ \ — J since Al-\-m is unity, 

"ill J b in 



and since s is unity, m = |a^, when m is estimated in the astronomical 
unit. 

* " Determination of the Seenlar Effects of Tidal Friction by a Grapliical Me- 
thod," " Proc. Roy. Soc," No. 197, 1879. 
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Again, since = ^Ma^, and since C is unifcj, therefore M=^/a^, where 
If is estimated in the mass unit. 

Therefore M'm/(M-\-m) is unity, when M and m are estimated in 
the mass unit, with the proposed units of length, time, and mass. 

According to the theory of elliptic motion, the moment of momentom 

of the orbital motion of the planet and satellite about their common 

Mm / 

centre of inertia is -jir-i — Oc^^l — e^. ISTow it has been shown that the 

factor involving Mis unity, and by (1) Qc^=:Q~i=:G^\ 

Hence, if we neglect the square of the eccentricity e, the m. of m. of 
orbital motion is numerically equal to Q ~ 3 or c^. 

Let i»=0~3=:c^. 

In this paper x, the moment of momentum of orbital motion, will 
be taken as the independent variable. In interpreting the figures 
given below it will be useful to remember that it is also equal to the 
square root of the mean distance. 

The moment of momentum of the planet's roiiation is equal to Cn ; 
and since is unity, n will be either the m. of m. of the planet's rotation, 
or the angular velocity of rotation itself. 

With the proposed units T=:f m/c^=:f a%~^, since w=f a^ ; and 
Cj;=|.^/a=flf/a^=fw. Mlma^=-f^pla. 

Also t^/q (a quantity which occurs below) is equal to ^^a^jvx^^. 

Now let t be the time, and let 2/ be the phase-retardation of the 
tide which I have elsewhere called the sidereal semi-diurnal tide of 
speed 2n, which tide is known in the British Association Ileport on 
Tides as the faster of the two K tides. 

Then if the planet be a fluid of small viscosity, the following are 
the differential equations which give the secular changes in the 
elements of the system : 

t=-i-«i"4/(l-"-) (2). 

dt 3 \ n / 

\ n J 



dt 


t Sill 

S 


di 
dt 


± -sm 



<-^){ 



l-^^^\ (4) 

n / 



_'??.= _iZ!8in4/ii±i) (5). 

dt g X 

1 de 1 "^^ • /I r 1 /n 18Q\ ^^x 

„. — .=. ;^— sin4jr.-( if — • I ..... {b). 

e dt g '^ x\ n J 

The first three of these equations are in eifect established in my 
paper on the " Precession of a Viscous Spheroid,"* § 17, p. 497, eq. (80). 

* " Phil. Trans.," Part II, 1879. 
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The suffix m^ to tlie symbols i and N there indicates that the equa- 
tions (80) only refer to the action of the moon, and as here we only 
have a single satellite, they are the complete equations. N is equal to 
ri/uQ, so that n^ disappears from the first and second of (80) ; also 
^=:1/s7^qQq^, and thus Uq disappears from the third equation. P = cos7*, 
Qrzrsin^', and, since we are treating i the obliquity as small, P=l, 
Qz=.i ; also \=^Q/n; the e of that paper is identical with the/ of the 
present one ; lastly f is equal to QqsQ-I, and since with our present 
units s-=-l, therefore /m d^/dt=^dQ~'^lnQdt^=^dxlnQdt. 

With regard to the transformation of the first of (80) into (4) of 
the present paper, I remark that treating i as small ^PQ—^^Q 
=z^i(l~2Qln), and introducing this transformation into the first of 
(80), equation (4) is obtained, except that i occurs in place of (i+j). 
Now in the paper on the '' Precession of a Viscous Spheroid " the in- 
clination of the orbit of the satellite to the plane of reference was 
treated as zero, and hence / was zero ; but I have proved in a paper 
^' On the Secular Changes in the Elements of the Orbit of a Satellite 
revolving about a tidally distorted Planet" (read before the Royal 
Society on December 18th, 1879, but as yet unpublished) that when 
we take into account the inclination of the orbit of the satellite, the P 
and Q on the right-hand sides of eq. (80) of '' Precession " must be 
taken as the cosine and sine of i-\-j instead of ^. Equations (5) and 
(6) are proved in § 10, Part II, and § 25, Part V of the unpublished 
paper, and the reader is requested to take them as established. 

The integrals of this system of equations will give the secular 
changes in the motion of the system under the influence of the 
frictional tides. The object of the present paper is to find an analy- 
tical expression for the solution, and to interpret that solution geome- 
trically. 

Prom equations (2) and (4) we have 



.dn. di 1 T^ . i /• 
^. — -j-n — = 4 — - sin4f 

dt dt ^^ -^ 



(^+i)+2ill ^ 



n 



But from (3) and (5) xd.jjdt-{-j dxjdt is equal to the same expres- 
sion ; hence 

>dn . di dj , . dx 

dt dt dt dt 

The integral of this equation is in-=:jx. 

'2' X , 1.^^ 

'J-'- ' ', • • • • • • e • \« /• 

Equation (7) may also be obtained by the principle of conservation 
of moment of momentum. The motion is referred to the invariable 
plane of the system, and however the planet and satellifce may interact 
on one another, the resultant m. of m. must remain constant in 
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direction and magnitude. Hence if we draw a parallelogram of 
which the diagonal is h (the resultant m. of m. of the system), and of 
which the sides are n and a?, inclined respectively to the diagonal at 
the angles i and j/, we see at once that 

sin?^ X 



smj n 

If i and J be treated as small this reduces to (7). 

Again the consideration of this parallelogram shows that 

'h^z=zn^-\-x^-\~ 2nx cos (?* -\-j) , 

which expresses the constancy of moment of momentum. If the 
squares and higher powers of i-\-j be neglected, this becomes 

}lz=:n-\-X (8). 

Equation (8) may also be obtained by observing that dnjdt-^- dxjdi'-^O, 
and therefore on integration n-\-x is constant. It is obvious from the 
principle of m. of m. that the planet's equator and the plane of the 
satellite's orbit have a common node on the invariable plane of the 
system. 

If we divide equations (4) and (6) by (3), we have the following 
results : — 

l^=L(l+iyiZL^^ (9). 

i dx 2u\ ij n — O 



(10). 



lde__l lln-18Q 

edx 2x n — Q 

But from (7) and (8) 

1+4=1+--—--, 

tAJ tAJ 

also 0=03"^, and n—h—'X. 

Hence (9) and (10) may be written 

d ^ . h 1 x^(h—x)—2 

I (~)p* 7 — _: t— 

dx 2 xQi—x) x^Qh—x) — ! 

d, 1 lW(h-x)-18 

— log e= — . — r--^ 

dx 2x x^(h—x) — l ^ 

h\x^(h-x)~2} li ,h x'- 

2x{h-x){x\h-~x)-l} x{h~x) 2x^-lix^ + l 

Therefore ■ 10^^=—+ + j • • • • \X^)' 

dx X h — x 2x'^'-]ix^ + l 



I 



(11). 



Also 



\\x^(li-x)-V 8 ^ 9 _ _ 7 x^{x-lv) 
2^H}i^x)—\) X 2' x^-lix^^\ 



260 



Mr. G. H. Darwin. 



[Mar. 18, 



Therefore 



d T 
- - log e- 

dx 



9 7 x'-(x~h) 



X 2 x'^ — JlQ^-^- I 



(la) 



These two equations are integrable as they stand, except as regards 
the last term in each of them. 

It was shown in a previous paper that the whole energy of the 
system, both kinetic and potential, was equal to -J- [t^^— -a?"^].* 

Then integrating (12) and (13), and writing down (7) and (8) 
again, and the expression for the energy, we have the following equa- 
tions, which give the variations of the elements of the system in terms 
of the square root of the satellite's distance, and independently of the 
time. 



X 



1 og i—\og —^ — 4- ^li 



h — x 



x^dx 



x^~hx^-{-l 



+ const. 



.^ ) 



loge— log.'c^— 1 



x^(x — h)dx , , 

— ^^ -/ — - + const. 



J 



# — hx^-^-l 



li — X. 

X 

n=:h — X. 



2E=^Qi—xy 



1 



X' 



J 



When the integration of these equations is completed, we si 
have the means of tracing the history of a fluid planet of small 
viscosity, attended by a single satellite, when the system is started 
Avith any given moment of momentum /i, and with any mean distance 
and (small) inclination and (small) eccentricity of the satellite's 
orbit, and (small) obliquity of the planet's equator. It may be re- 
marked that li is to be taken as essentially positive, because the sign 
of }i merely depends on the convention which we choose to adopt as 
to positive and negative rotations. 

These equations do not involve the time, but it will be sbown later 
how tbe time may be also found as a function oi x. It is not, how- 
ever J necessary to find the expression for the time in order to know 
the sequence of events in the history of the system. 

Since the fluid which forms the planet is subject to friction, there- 
fore the system is non- conservative of energy, and therefore x must 
change in such a way that E may diminish. 

If the expression for E be illustrated by a curve in which 'E is the 
vertical ordinate and x the horizontal abscissa, then any point on this 
" curve of energy " may be taken to represent one configuration of 
the system, as far as regards the mean distance of the satellite. 



* '' The Secular Eteocts of Tidal Friction," 
1879, eq. (4). 



•5 



^'Proc. Eoy Soc.," Fc. 19 
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Tlien suck a point; must always slide down a slope of energy, and we 
shall see which way x must vary for any given configuration. This 
consideration will enable us to determine the sequence of eventSj 
when we come to consider the expressions for i^ e, j^ n in terms of x. 

We have now to consider the further steps towards the complete 
solution of the problem. 

The only difficulty remaining is the integration of the two expressions 
in the first and second of (14). From the forms of the expressions 
to be integrated, it is clear that they must be split up into partial 
fractions. The forms which these fractions will assu'me will of course 
depend on the nature of the roots of the equation x^~-lix^-\-\-=:-0. 

Some of the properties of this biquadratic were discussed in a 
previous paper, but it will now be necessary to consider the subject 
in more detail. 

It will be found by Ferrari's method that 

where X^— 4X.~/i^=0. 

By using the property (X'^-^li) (\t + /0^^^? ^^^^ expression may be 

written in the form 

[.r;^4.i(Xt~/j.)}3-f {|(xl-^) v/r+2M^}3J X 

Which is of course equivalent to finding all the roots of the bi- 
quadratic in terms of h and \. 

ETow let a curve be drawn of which lir' is the ordinate (negative 

Fig. 1. 









1 








1 



The ordinates are dmwn to one-third of the scale to which t3ie abscissa) are drawn. 

VOL. XXX. T 
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values of hP' being admissible) and X the abscissa ; it is shown in fig. 1. 
Its equation is /?'^-.=X(X^— 4). 

It is obvious that 0A=0A'=2. 

The maximnni and mininiiiin values of JiP (viz., Bh, B'b') are given 
•bjS\^=4^oT\=±2iSK 

Then B5r=B^&'=::-23/3t + 4. 2/3i=:(4/3^)3. 

Since in the cubic, on which the solution of the biquadratic depends, 
h^ is necessarily positive, it follows that if Ji be greater than 4/3' the 
cubic has one real positive root greater than OM, and if h be less than 
4/3% it has two real negative roots lying between and OA', and 
one real positive root lying between OA and OM. 

To find OM we observe that since li^ is equal fco (4/3^)^, and since 
the root of X^— 4X — h^=-0 which is equal to — 2/3^ is repeated twice, 
therefore, if e be the third root (or OM) we must have 

X + |^V--e)^X3-4X™^^); 

whence (2/3^)3err: (4/3^2)2^ and e or OM=4/3l 

iSTow 0A=:2 ; hence, if h be less than 4/3^, the cubic has a positive 
root between 2 and 4/3% and if li be greater than 4/3', the cubic has a 
positive root between 4/3- and infinity. 

It will only be necessary to consider the positive root of the cubic. 

Now suppose h to be greater than 4/3^ 

Then it has just been shown that X is greater than 4/3% and hence 
(X being positive) 3X^ is greater than 16X, or 4(X.^~4X) greater tlian 
X^, or 47z^ greater than X*^, or 2/i\'"l greater than unity. 

Therefore {|()a4-/t) v^l^^X^P^^ ---{K^+A) v^m"l-~lp. 
Thus the biquadratic has two real roots, which we may call a and 1). 
Then a=:|(Xl + 7?.)[l4- ^^mx^h~ll 

b =r K^^ + ^0 [1 - ^2lx=f^] . 
It will now be proved that a is greater and b less than |/i. 
l^ow a> or <|/% 

as (Xl + 7i)[l+v^2^i=l]> or <3/% 

(XI H-^) / 

as ^^^27t--Xl> or <2^-XI, 

as XI^-7^> or <X^v^2F^, 
as XS + 2^.XI+7z.2> or <2/at--X'\ 
as 2\^-i-h^> or <0. 
Since the left hand side is essentially positive, a is greater than f /?.. 
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Again 1)> or <~^h, 

as (U + li)\l—'^2}h)F¥^r]> or <3^, 

as -~~|-^V2/^-Xt> or <2^--XI. 

Since tlie left-hand side is negative and the right positive, the left is 
less than the right, and therefore b is less than ^Ji. 
If, therefore, Ji be greater than 4/3^, we may write 

where a— f ^, f ^— b are positive, and where a is negative. 

We now turn to the other case and suppose h less than 4/8*. All 
the roots of the biquadratic are now imaginary, and we may put 

If a be taken as — ^(Xt— ^), then ey is }(Xt + ^). 

Then it only remains to prove that 7 is greater than f ^. 

Now 7 > or < ffh, 

as Xt > or < 2]i, 

asX''^> or <4P=4(X^— 4X), 

as 16 > or < 3X^, 

as 4/3^ > or < X, 

but it has been already shown that in this case, X is less than 4/3\ 
wherefore 7 is greater than ^h. 



We may now proceed to the required integrations. 
First case where h is greater than 4/3^ 

Let a34--/^a3Hl = (^™a)(a;--b)[(a?-a)3 + /33], 

so that the roots are a, b, oc+/3i. 

Also let a be the root which is greater than f /t, b that which is less^ 
and let 



/jjArf 



To find the expression for i we have to integrate -- 



a?* — /?.aj^4-l 

Let /(a?) = («— a)-^(0), and let a5^//(a3)=J./(a?— a) +(p(w)l^{m) 
Then a?^ (x — a) = Af(x) + (x — a) ^^ (a?) . 

Hence J_=a^//'(^)- 

If, therefore, f(x):=x'^—]ix^-{-l, J.=l/(4a— 3/^)=l/4a|. 

T 2 



% 



.ij OS- 
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Thas the partial fractions corresponding to tlie roots a and b are 



1 



4a|aj — a ^b^a^ — b 



(15). 



If tlie pair effractions corresponding to the roots oc^(3i be formed 
and added together, we find 

1 — x^ix — a)-f/3^ 



2(«i2 4-/33) [(ic-^) 3-^/33] 



The siiiii of (1,5) and (16) is eqnal to 



X'' 



x^ — hx^ + 1 



5 and 



9'^dx 



x^~li9?-^\ 4a: 



log {X00 2b)— l0g(a3 CO b) — ^^ 



+ 



4bi 



— arc tan 
2(«i3-f/33) 8 



4(«i3^/33 

X — oc 



log[(^-^)3-f/33J 



(17), 



Substituting in the first of (14) we have 



%^=^. 



X 



(xoosb)8nx exp. 



hp 



X — cc 



- arc tan 

_4(^i^ + /3'^) /3 



h — ''A 



h ha .^ 

(.TC0b)81^ [(OJ— <%)''^-}-/53]8{ai=T3^ 



(18), 



where ^ is a constant to be determined by the value of ?', which 
corresponds with a particular value of x. 

From the third of (14) we see that by omitting tbe factor 
xl(Ji~-93) from the above, we obtain the expression for y. 

To find the expression for e we have to integrate -L-A^ Jl , 

x'^ — hx^ + l 



x^(x-h)=l(4^x^~-Slix'-) 



17 '> 



and therefore 

x^(^d — Ji)dx 



X^-Jliff^ + l 



^log (x'^—]ix^-\-l)—^h 



X 



^dx 



xJ^ —lix^ -{-X 



The integral remaining on the right hand has been already deter 
mined in (17), Then substituting in the second of (14), we have 



Q^s 



Bx 



e: 



(»;'i-/^;»Hiy 



h 



(a3coa)8ai exp. 



A/3 



4(«l3 4_^3) 



arc tan 



X—cx, 



-^\ 



^ 



V. 



h fto^i 

((»(X)b)8bi [(a? — a)^- + ^3-|8(a-/-'+i3^ 



(19),- 



/ 



where B is a constant to be determined by the value of e, correspond 
ing to some particular value of x. 
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From this equation we get the curious relationship 



265 



e- 



B 



X' 



9 



Ai' (m'^-hx^ + l)^ 



f 



fi 



20). 



This last result will obviously be equally true even if all the roots 
of #—^#+1=0 are imaginary. 

In the present case the complete solution of the problem is com- 
prised in the following equations : — 



(r(j(X)a)8aiexp. 



J 



■=A: 






\ 



h ha t 



X 



e 



h — X 



J' 



> 



# 



■nzzzh—x. 



2Bz=z(h-xy~-±^, 



•7 






It is obvious that the system can never degrade in such a way 
that m should pass through one of the roots of the biquadratic 
.^;i._/ia;3^][=0. Hence the solution is divided into three fields, viz., 
(i) a?=+oo to a?=a; here we must write a?— a, a?— b for the cecoa, 
c^cob in the above solution; (ii) a?=atoa?=b; here we must write 
a-aj, a?— b (this is the part which has most interest in application to 
actual planets and satellites) ; (iii) a5=:btoa?= — oo; here we must 
write a— a?, h—x. When x is negative the physical meaning is that 
the revolution of the satellite is adverse to the planet's rotation. 

By referring to (4) and (6), we see that i must be a maximum or 
minimum when w=20, and e a maximum or minimum when *^= |f Q. 
Hence the corresponding values of x are the roots of the equations 
tij*— /iaj3 + 2=0, and x^—lhx^-\-\l^O respectively. 

Since 






1 



1 



•it; 



■o!,^(x — a) -\- ji'^ 



x^-hd^-\-l 4aire-a 4b^x~h 'Z{oc{-+fi^) [(aj-^)3 + /33] 

Therefore 

1 



1 



4a, 



,(,^™|>)[(,,-.^)34-^3]^^(a'.-a)[0^-a)2 + ^^] 






ifcD| 



2(ci5^+j83) 



[— iKi(a3— a) +/3^](a}~-a)0(;-~- b). 
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Hence the coefficient of x^ on the right-hand side must be zero, and 

therefore — — tct— tt^^ — -_=0, 

4ai 4bi 2(«i3 + ^3) 



And __— =—--}- 



8ai 8bi 4(<%i2 + ^2) 



ISTow when q3=:+ oo, arc tan -^-^ = i^3 and when 33=:=— oo, it is 

[3 

equal to — Jtt. 

Hence when 03=:+ go, j;'=: J. exp. [i'7"/^/5/8(a|^H-/33)], i^z—j^^ the 

upper sign being taken for + 00 and the lower for — oo. 

Then since j tends to become constant when 03=+ ^? and since 
9 — 1= ^J^, therefore when «? is very large e tends to vary asa3V. 

If X be very small j has a finite value, and i varies as x, and e varies 

j, i, and e all become infinite when a^— b, and x also becomes infinite 
when x=:Ji. 

This analytical solution is so complex that it is not easy to under- 
stand its physical meaning ; a geometrical illustration will, however, 
make it intelligible. 

The method adopted for this end is to draw a series of curves, the 
points on which have x as abscissa and i, j, e, n, ]]] as ordinates. The 
figure would hardly be intelligible if all the curves were drawn at once, 
and therefore a separate figure is drawn for i, j, and e ; but in each 
figure the straight line which represents qi is drawn, and the energy 
curve is also introduced in order to determine which way the figure is 
to be read. The zero of energy is of course arbitrary, and therefore 
the origin of the energy curve is in each case shifted along the vertical 
axis, in such a way that the energy curve may clash as little as 
possible with the others. 

It is not very easy to select a value of Ji which, shall be suitable for 
drawing these curves within a moderate compass, but after some con- 
sideration I chose h=.2'6, and figs. 2, 3, and 4 are drawn to illustrate 
this value of h. If the cubic X-^—4\ — (2- 6)^=0, be solved by Cardan's 
method, it will be found that X=2'b741, and using this value in the 
formula for the roots of the biquadratic we have 

a34_2-6arHl = (a3-2-539)(a3--826)[(a3 + -382)3+ (-.575)3]. 

Hence a = 2-539, b==-826, ^:=--382, I3=:'576, |/^=::l-95, and 
4.a3=:2-356, 4^\=:4^'4^96, «i =.2-332, ocY^-h(3^=b'771. 
Then we have 
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.__ , (2-539 cvD aj)"S53 exp. ['062 arc tan (1740a)+ '665)] ^ 
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The maximum, and minimum values of i are given by the roots of 
the equationa?*— 2'6a)^H-2=:0, viz., a? =2*46 7 and !e=l'103. The max- 
imum and minimum values of e are given by the roots of the equation 
;^4 _ 2*603 +^|-=0, vijs., a?=:2*495 and a)=l-0095. The horizontal asymp- 
fcotes for i\A and ^\A are at distances from tlie axis of x equal to 
exp. (•062x^7r) and exp. ( — •062x^7r), which are equal to 1'102 and 
•908 respectively. 

Fig. 2 shows the curve illustrating the changes of ^, the obliquity of 
the equator to the invariable plane. 

The asymptotes are indicated by broken lines ; that at A is given 
by 03= '826, and is the ordinate of maximum energy ; that at B is given 
by a) =2*6, and gives the configuration of the system for which the 
planet has no rotation. The point C is given by a? =2*5395 ^"^^ lies on 
the ordinate of minimum energy. Geometrically the curve is divided 
into three parts by the vertical asymptotes, but it is further divided 
physically. 

The curve of energy has four slopes, and since the energy must 
degrade, there are four methods in which the system may change, 
according to the way in which it was started. The arrows marked on 
the curve of obliquity show the direction in which the curve must be 
read. 

Since none of these four methods can ever pass into another, this 
figure really contains four figures, and the following parts of the figure 
are quite independent of one another, viz. : (i) from ~ oo to O ; (ii) 
from A to O ; (iii) from A to C ; (iv) from + oo to C. The figures 
3 and 4 are similarly in reality four figures combined. For each of 
these parts the constant A. must be chosen with appropriate sign ; but 
in order to permit the curves in fig. 2 to be geometrically continuous 
the obliquity is allowed to change sign. 

The actual numerical interpretation of this figure depends on the 
value of A. Thus if for any value of x in any of the four fields the 
obliquity has an assigned value, then the ordinate corresponding to 
that value of x will give a scale of obliquity from which all the other 
ordinates within that field may be estimated. 
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Diagram for Obliquity of Planet's Equator. — First case. 

As a special example of tMs we see that, if tlie obliquity bo zero at 
any point, a consideration of the curve will determine wbether zero 
obliquity be dynamically stable or not ; for if the arrows on the curve 
of obliquity be approaching the axis of x, zero obliquity is dynamically 
stable, and if receding from the axis of x, dynamically unstable. 

Hence from a3=4-oo to B^ zero obliquity is dynamically unstable, 
from -— 00 to O and A to O dynamically stable, and from A to B, first 
stable, then unstable, and finally stable. 

The infinite value of the obliquity at the point B has a peculiar 
significance, for at B the planet has no rotation, and being thus free 
from what Sir William Thomson calls '' gyroscopic domination," the 
obliquity changes with infinite ease. In fact at B the term equator 
loses its meaning. The infinite value at A has a different meaning. 
The configuration A is one of maximum energy and of dynamical 
equilibrium, but is unstable as regards mean distance and planetary 
rotation ; at this point the system changes infinitely slowly as regards 
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time, and therefore tke infinite valne of tlie obliquity does not indicate 
an infinite rate of ctange of obliquity. In fact if we put .*=Q in (1) 
we see that c^^/iJc^if =--^(7^/2) sin 4/. However, to consider this case 
adequately we should have to take into account the obliquity in the 
equations for (Injdt and dcc/dt, because the principal semi-diurnal tide 
vanishes when %=Q. 

Similarly at the minimum of energy the system changes infinitely 
slowly, and thus the obliquity would take an infinite time to vanish. 

We may now state the physical meaning of fig. 2, and this interpre- 
tation may be compared with a similar interpretation in the paper on 
" The secular effects of tidal friction," above referred to. 

A fluid planet of small viscosity is attended by a single satellite, and 
the system is started with an amount of positive moment of momentum 
which is greater than 4/3^, with our present units of length, mass and 
time. 

The part of the figure on the negative side of the origin indicates a 
negative revolution of the satellite and a positive rotation of the 
planet, but the m. of m. of planetary rotation is greater (by an amount 
h) than the m» of m. of orbital motion. Then the satellite approaches 
the planet and ultimately falls into it, and the obliquity always 
diminishes slowly. The part from O to A indicates positive rotation 
of both parts of the system, but the satellite is very close to the planet 
and revolves round the planet quicker than the planet rotates, as in 
the case of the inner satellite of Mars, Here again the satellite 
approaches and ultimately falls in, and the obliquity always diminishes. 

The part from A to indicates positive rotation of both parts, but 
the satellite revolves slower than the planet rotates. This is the case 
which has most interest for application to the solar system. The 
satellite recedes from the planet, and the system ceases its changes 
when the satellite and planet revolve slowly as parts of a rigid body — 
that is to say, when the energy is a minimum. The obliquity first 
decreases, then increases to a maximum, and ultimately decreases to 
zero.* 

The parti from infinity to C indicates a positive revolution of the 
satellite, and from infinity to B a negative rotation of the planet, but 
from B to a positive rotation of the planet, which is slower than the 
revolution of the satellite. In either of these cases the satellite 
approaches the planet, but the changes cease when the satellite and 
planet move slowly round as parts of a rigid body — ^that is to say, when 
the energy is a minimum. If the rotation of the planet be positive, 
the obliquity diminishes, if negative it increases. If the rotation of 
the planet be nil, tbe term obliquity ceases to have any meaning, since 
there is no longer an equator. 

* According to the present theory, the moon, considered as being attended by the 
vsarth as a satellite, has gone through these changes. 
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Fig. 3 illiistrates tlie changes of inclination of the satellite's orbit, 
and may he interpreted in the same way as &g. 2. It appears from the 
part of the figure for Avhich w is negative, that if the revolution of the 
satellite he negative, and the rotation of the planet positive, but the 



Fig. 3, 




Diagram for Inclination of Satellite's Orbit. — First c?ise. 



m. of m.. of planetary rotation greater than that of orbital motion^ 
thenj as the satellite approaches the planet, the inclination of the orbit 
increases, or zero inclination is dynamically nnstable. In every other 
case the inclination will decrease, or zero inclination is dynamically 
stable. 

This result undergoes an important modification Nvhen a second 
satellite is introduced, as will appear in the unpublished paper. 

Mg. 4 shows a similar curve for the eccentricity of the orbit. The- 
variations of the eccentricity are very much larger than than those of 
the obliquity and inclination, so that it was here necessary to draw the 
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orcliiiates on a much reduced scale. It was not possible to extend the 
figure far in either direction, because for large Talues of «, e varies as 
a high power of x (viz., V")* The curve presents a resemblance to that 
of obliquity, for in the field comprised between the two roots of the 
biquadratic (viz., between A and 0) the eccentricity diminishes to a 
minimum, increases to a maximum, and ultimately vanishes at 0. 



Fig. 4. 




Diagram for Eccentricitj of Satellite's Orbit. — First case. 



This field represents a positive rotation both of the planet and satel- 
lite, but the satellite revolves slower than the planet rotates. This part 
represents the degradation of the system, from the configuration of 
maximum energy to that of minimum energy, and the satellite recedes 
from the planet, until the two move round slowly like the parts of a 
rigid body. 

In every other case the eccentricity degrades rapidly, whilst the 
satellite approaches the planet. 

The very rapid rate of variation of the eccentricity, compared with 
that of the obliquity would lead one to expect that the eccentricity of 
the orbit of a satellite should become very large in the course of its 
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evolution, whilst the obliquity should not increase to any very large 
-extent. But it must he remembered that we are here only treating a 
planet of small viscosity, and it will appear, in the unpublished paper 
above referred to, that the rate of increase or diminution of the eccen- 
tricity is very much less rapid (per unit increase of x) if the viscosity 
be not small, whilst the rate of increase or diminution of obliquity 
(per unit increase of x) is slightly increased with increase of viscosity. 
Thus the observed eccentricities of the orbits of satellites and of 
obliquities of their planets cannot be said to agree in amount with the 
theory that the planets were primitively fluids of small viscosity, 
though I believe they do agree with the theory that the planets were 
fluids or quasi-solids of large viscosity. 



We now come to the second case, where li is less than 4/3*. The 
biquadratic having no real roots, we may put 

It has already been shown that a is negative, and 7 greater than ^h. 

Let a = ^h — ai, fy = ry^-|-|j^. 

Then by inspection of the integral in the first case we see that 



exp. -~~ — J- -^^ arc tan 4- — - — -- — ™ arc tan — ^ 

l4(£Kj'"-|-^^) ji '^('Yi "f"^ ) ^^ 

The rest of equations (21), which express the other elements in 
iterms of j and ro, remain the same as before. 
By comparison with the first case, we see that 

x^ __ 1 — «j(re— a) +/3^ 1 72(^3—7) 4- cl3 



x'^—hx^ + l 2(^3 '■+/i^) i^Jc.—ocY+fi^ 2(7/ + f33) (a?— 7)34-^2 

•On multiplying both sides of this identity by f^* — /^a?'M-lj and equating 
the coefficients of rc'^, we find 



0= ~"i - + "/I 



Therefore 






'Thus when x is e€|iial to +00 



J =^4 exp. ± 



8(^i3^^2) 8(7/4-^3). 
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the upper sign being taken for +00, and the lower for —00. This 
expression gives the horizontal asymptotes for j and i. 

In order to illustrate this solution, I chose h^^^l, and found by tri- 
sronometrical solution of the cubic X'^— 4X — 1=0, X:=:r 2-1149, and 
thence 
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Diagram for Ohliqiiity of Planet's Equator. — Second case. 



Tliese solutions are illustrated as in the previous case by the three 
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figures 5, 6, 7. There are here only two slopes of energy, and hence 
these figures each of them only contain two separate figures. 

Fig. 5 illustrates the changes of ?*, the obliquity of the equator to 
the invariable plane. 

In this figure there is only one vertical asymptote, viz., that corre- 
sponding to a?=l. For this value of x the planet has no rotation, is 
free from "gyroscopic domination," and the term equator loses its 



meaning. 



The figure shows that if the rotation of the planet be negative, but 
the m. of m. of planetary rotation less than that of orbital motion, 
then the obliquity increases, whilst the satellite approaches the planet. 

This increase of obliquity only continues so long as the rotation of 
the planet is negative. The rotation becomes positive after a time, and 
the obliquity then diminishes, whilst the satellite falls into the planet. 
In the corresponding part of iig. 2 the satellite did not fall into the 
planet, but the two finally moved slowly round together as the parts of 
a rigid body. 

If the revolution of the satellite be negative, and the rotation of 
the planet positive, but the m. of m. of rotation greater than that of 
revolution, the obliquity always diminishes as the satellite falls in to 
the planet. 



Fia. 6. 




Diagram for Inclination of Satellite's Orbit. — Second case. 
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Figs. 2 and 5 only differ in tlie fact tliat in the one there is a true 
maximum and a true minimum of obliquity and energy, and in the 
other there is not so. In fact, if we annihilate the part between the 
vertical asymptotes of fig. 2 we get fig. 5. 

Fig. 6 illustrates the changes of inclination of the orbit. It does 
not possess very much interest, since it simply shows that however the 
system be started with positive revolution of the satellite, whether the 
rotation of the planet be positive or not, the inclination of the orbit 
slightly diminishes as the satellite falls in. 

And however the system be started with, negative revolution of the 
satellite, and therefore necessarily positive rotation of the planet, the 
inclination of the orbit slightly increases. Fig. 6 again corresponds 
to fig. 3, if in the latter the part lying between the maximuin and 
minimum of energy be annihilated. 



Fig. 7. 




Diagram for Eccentricity of Satellite's Orbit. — Secoud case. 

Fig. 7 illustrates the changes of eccentricity, and shows that ifc 
always diminishes rapidly however the system is started, as the satellite 
falls into the planet. This figure again corresponds with fig. 4, if in 
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tbe latter -i-lie parts between the niaximiini and mininmm of energy be 
annihilated. 

These three figures may be interpreted as giving the varions stabilities 
and instabilities of the system, just as was done in the first case. 



The solution of the problem, which has been given and discussed, 
above, gives merely the sequence of events, and does not show the rate 
at which the changes in the system, take place. It will now he shown 
how the time may be found as a function of x. 

Consider the equation 

/"//>' T'-' I O \ 

— =-i- — sm 4/1 1— — . I, 
at g \ n / 

/is here the angle of lag of the sidereal semi-diurnal tide of speed 2?2. ; 
then by the theory of the tides of a viscous spheroid, tan2/=2??/p, 
where p is a certain function of the radius of the planet and its density, 
and which Taries inversely as the coefficient of viscosity of the 
spheroid.* 

Since by hypothesis the viscosity is small, / is a small angle, so that 
sin 4if may be taken as equal to 2 tan 2/. Thus, sin 4//?2. is a constant, 
depending on the dimensions, density, and viscosity of the planet. 

It h{:is already been shown that t^ varies as 03"^^, and g is a constant, 
which depends only on the density of the planet. Hence, the above 
equation may be written 

where K is a certain constant, which it is immaterial at present to 
evaluate precisely. 

Since n^i—9} and 0=aj"'\ we have 

Kdt: ^ ^^ 



or K/=: -— I — ^ — -f a const. 

Jx^—hx-^ + l 

The determination of this integral presents no difficulty, but the 
analytical expression for the result is very long, and it does not at 
present seem worth while to give the result. The actual scale of time 
in years will depend on the value of K, and this is a subject of no 
interest at present. 

It will, however, be possible to give an idea of the rate of change 
of the system without actually performing the integration. This may 

* "On the Bodily Tides of Yiscous and semi-elastic Spheroids," &c. "Phil. 
Trans.,," Part I, 1879, p. 13, § 5. 
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be done by drawing a curve in wbicli the ordinates are proportional 
to dt/dx, and the abscissae are x. The equation to this curve is then 






dt 



■X 



15 



div # — ha^-i"! 

The maximum and minimum values (if any) of cU/d^^ are given by the 
real roots of the equation 

One of such roots will be found to be intermediate between a and b, 
and the other greater than a. 
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Diagram illustrating the Rate of Change of the System. 

Fig. 8 shows the nature of the curve when drawn with the free 
hand. It was not found possible to draw this figure to scale, because 
when ^=2'6 it was found that the minimum M was equal to "SS, and 
could not be made distinguishable from a point on the asymptote A, 
whilst the minimum m was equal to about 900,000, and could not be 
made distinguishable from a point on the asymptote C. 

The area intercepted between this curve, the axis of x, and any pair 
of ordinates corresponding to two values of a?, will be proportional to 
the time required to pass from the one configuration to the other. 

Where dl/dx is negative, that is to say, when the satellite is falling 
into the planet, the areas fall below the axis of x. This is clearly 
necessary in order to have geometrical continuity in the curve. 

The figure shows that the rate of alteration in the system becomes 

VOL. XXX. U 
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very slow wlien tlie satellite is far from the planet ; this must indeed 
obviously be the case, because the tidal effects vary as the inverse sixth 
power of the satellite's mean distance. 



V. *^ On the Modifications of the Spectrum of Potassium which 
are Effected by the Presence of Phosphoric Acid, and oxi 
the Inorganic Bases and Salts v^hich are found in combina- 
tion with Educts of the Brain." By J. L. W. Thudighum, 
M.l)., F.R.O.P.L. Communicated by John Simon, C.B.. 
F.R.S., &c. Received March 10, 1880. 

Among the results of a large investigation on which I have for 
many years been engaged in regard of the chemistry of the brain, I 
had been led to conclude that the so-called "protagon" of Oscar 
Liebreich is not a definite chemical body, but is a variable mixture of 
several bodies. This conclusion of mine (which agrees with opinions 
expressed on the same subject by Strecker, Diaconow, and Hoppe- 
Seyler) was published by me in 1874,* and endeavours to controvert 
it have since then been made, on several occasions, by Dr. Arthur 
Gamgee.t; Last summer, he brought before the Royal Society J his 
contentions for the chemical individuality of " protagon "; and it 
fortunately was in my power shortly afterwards to publish evidence, 
which, I believe, those who will take the trouble to follow it will 
find quite unanswerable, that Dr. Gamgee's contentions were mis- 
taken. § Part of my evidence to that effect consisted in showing by 
quantitative analyses that Dr. Gamgee's so-called ^'protas'on" con- 
tains 0*7 per cent, of potassium; secondly, that in connexion with 
trifling differences in the extraction process, the proportion of potas- 
sium in different specimens of " protagon " can be made to range 
from a trace to 1*6 per cent. ; thirdly, that with the variable quantities 
of potassium the quantities of phosphorus and other ingredients will 
also vary. 

In the last published number, No. 200, p. Ill, of the " Proceedings 
of the Royal Society," I find that Dr. Gamgee has recently brought 
the question again under notice of the Society, and that, in doing so, 
he especially rests his case upon the following statement made by his 
colleague, Professor Roscoe, on the subject of some examinations, 

* " Reports of the Medical Officer of the Privy Council and Local Grovernment 
Board." New Series, No. III. 

+ " Zeitschrift £ur Physiol, Chemie," vol. iii, p. 260 ; ^' Ber. Dentsch. Chem. 
aes.," 1879, &c. 

X ''Proc. Roy. Soc," vol. xxix, p. 151. 

§ " Annals of Chemical Medicine." Edited by J. L. W. ThndichuiD. Yoi. i, 
p. 254 
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Diagram for Efcentricitv of Sutellilc':i Orbit. — Fii-st cast-. 
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Diagrum for Obliquity of Planet's Equator. — St-cond case. 
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DiayraiH for Jneliiiatioii of iSak'lHte's Orbit.— Second case. 
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DiayraiH for Jneliiiatioii of iSak'lHte's Orbit.— Second case. 




Diagram illustrating the Kate of Change of the System. 



